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ABSTRACT 


Integration of photovoltaic systems (PVs) without battery energy storage into distributed grids is considered in this paper to 
evaluate of this generation to their operating parameters. System structure (including power converter, measurement and control 
system) and working modes of PVs will be introduced and analyzed to show the its power variation. The system of bus power 
balance and Newton-Raphson algorithm is also established in this paper to determine operating parameters of the microgrid. A 35- 
kV distributed grid is simulated by the ETAP software corresponding to a fixed value for load buses and varied values of power of 
solar irradiance. Simulation results showed the change of power flows in branches and the increase of bus voltage in cases of 
having power participation of PVs at a bus (enough power of solar irradiance). Due to working as a distributed generator, research 
results are useful for dispatchers to make operating plans for whole power system and implement PVs at other buses or change 


transmission lines with suitable sections. 
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INTRODUCTION 


PVs is a system harnessing power from photovoltaic generator (PVG) through power converters. The integration of PVs into 
distributed grids can affect to operating system parameters. It must be evaluated by using mathematical and simulation tools to 
make suitable operating plans for whole power system at any time [1]. 

Power participation of PVs depends on the capability of photovoltaic conversion corresponding to each type of PVG, power of 
solar irradiance (G), temperature (T) at p-n junction and control method [2]. To harness power from PVs and connect to the power 
system, power converters play an important role to regulate PVG [3]. They can help to harness maximum power from PVG, 
integrate PVG and energy storage (option), synchronize to the grid and generate power to the grid. At each time, dispatchers must 
find out characteristics of PVs to exactly determine integrating power between PVs and the grid as nature or requirement that help 
to hold stability for whole system. 

To have information about power flows in branches and bus voltage, Newton-Raphson method was proposed to solve the 
system of power balance equations in power system. Approximate solutions can be determined by classifying type of buses and 
using impedance matrix. Because the value of solutions is near to rated value, this method has a fast speed of convergence and is 
often chosen to determine operating parameters of power system by computer [4], [5]. 

Corresponding to above analysis, it must be made clear about the method to harness power, the integration between PVs and 
the grid at any time and the effect of PVs into operating parameters of distributed grids. To execute this purpose, the next section 
will represent the structure and characteristic of system harnessing PVG. Section III will introduce contents of Newton-Raphson 
method applied to determine operating parameters of the grid. Section IV will show some simulation results for a grid example in 


some scenarios of G. The last section will give out some conclusions. 


POWER STRUCTURE AND CONTROL FUNCTIONS OF PVs 


In grid-connected systems, PVG can be harnessed and synchronized by DC/DC and DC/AC converters. Power structure and control 
functions of PVs are described in Fig. 1 [3]. 


PV power PV =) er PV =) er PV power 


= Cc) 


Filter 
Pe Siee 


Current/Vo a ge e4 Grid = onization 
contr a contr e4 
MPPT Harmonic [eeeeee ae | Power 
Compensation [eeeeee ae | Control 


> 
a 


v . . 
<i(__)e Solar irradiance 
& Transformer 


Li 


Control functions 


Fig. 1. Power structure and control functions of PVs 


Power structure of this system has some main units such as Cac capacitor, DC/DC and DC/AC power converters, filter and 
transformer. Cac capacitor helps to charge and discharge energy instantaneously that contributes to hold stable voltage at DC bus. 
The DC/DC converter regulates PVG to harness maximum power from this generation at any time. The DC/AC helps to 
synchronize with the grid and generate PV power to the grid. The filter and transformer play an important role to filt high order 
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harmonic of the current integrating between PVs and the grid. Output voltage of the transformer is also suitable for the requirement 


of the grid. Due to having no energy storage and small power loss, almost PV power runs through the DC/DC and DC/AC to the 
grid without charging. For this reason, power from PVs always fluctuates corresponding to the variation of G. 

Some main measurement quantities are voltage and current at output terminals of PVG; voltage at DC bus; voltage and current 
at point of common coupling with the grid. These quantities must be collected in this structure to execute control functions, such as 
control current or voltage at output terminals of PVG to operate it at maximum power point; regulate voltage at DC bus at a 
reference value; control to compensate high order harmonic or regulate power integration between PVs and the grid. 

Although power generating from PVG depends on many factors but it only varies much in case of much variation of G. Peak 
points, called maximum power point (MPP), moves up or down instantaneously and immediately corresponding to the change of G 
as described in Fig. 2 [2]. 


Power 


Fig. 2. Movement of MPP corresponding to the change of G 


Power integrating between PVs and the grid is zero or near zero when G is very small. It also reached to rated value at any time in 
daytime if standard test condition exists. Due to always fluctuation and large range variation, power generating from PVs can affect 
much to the grid and must be evaluated more detailed through simulation process before having decisions to dispatch whole 
system. 


EFFECT PVs WITHOUT ENERGY STORAGE TO POWER FLOWS AND BUS VOLTAGE IN 


DISTRIBUTED GRIDS 
A simple grid with the participation of PVs at bus B is considered in Fig. 3. 
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Fig. 3. A simple power system 


In this system, active and reactive power losses in transmission line can be determined by (1) [5]: 
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where, Pt and APiine (kW) are active power of load and active power loss in transmission line, 


Ppv (kW) is the integration power between PVs and the grid, 
Qt and AQtine (kVAr) are reactive power of load and reactive power loss in transmission line, 
Ua (kV) is voltage at bus A, 


Riine and Xiine (Q) are resistance and reactance of the transmision line. 


Voltage at bus B can be determined by (2): 
= (P, ~ Po JR tine + Q: Xiine 


AU, 
line U, 


©00 


At night-time (without solar irradiance), power integrating between PVs and the grid is zero. This case is same as cases without 
implementing PVs at bus B. At daytime (having solar irradiance), PVs will help to reduce power mobilized from source, power and 
voltage losses (APpv and AUiine). 

When the grid becomes more complex with many branches and buses, power flows in branches and bus voltage are also more 
difficult than above simple grid. In these case, Newton-Raphson method can be used to solve this problem due to fast convergence 
and high accuracy [4]. It means that the active effect of PVs to the grid will be made easier and clearer if it is implemented at many 
buses with suitable capacity. 

To determine operating parameters for N-bus grid by using Newton-Raphson method, system of power balance equations at 
the i bus can be defined by (3) and (4) [5]: 


N 

U’y,, cosw, + YUU, cos(8,-5,-y,,)=AP, eee 
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where: i1=1,N; AV, = Pe =Fas AQ, = Q5 Oe U; = U73;; Yj a YiZVj 
Pui and Qui are active and reactive load power at the i bus; Pci and Qci are active and reactive power of PVs at the i bus. 


Because PVs is considered as an active generator, integrating power between PVs and the grid at any time depends on G and T as 
defined by (5): 


©OO 


In this problem, solutions are module and angle of bus voltage at the i bus. Solutions at the (k+1) step can be determined by 


(6) [5]: 
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where, AUi® and Adi are module and angle errors of bus voltage at the the i bus and the k"* step. These errors can be calculated 
by (7) at the k® step [5]: 


(k) (k) 
As] __[ AP 
= OOO 


Aue, AQ” 


where, J! is the inversed matrix of Jacobian matrix, 
OPi®)/06i®, EPI /OUi, CQi®)/06i®) and OQ: /OUi® are elements of Jacobian matrix. 


APi® and AQi® are active and reactive power errors at the k" step 


Newton-Raphson algorithm applied to determine operating parameters with the participation of PVs is described in Fig. 4 [5]. 
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Fig. 4. Newton-Raphson algorithm to determine operating parameters in the grid 
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SIMULATION RESULTS 


The diagram of the example grid with the participation of PVs at bus P4 is considered in Fig. 5. 


LoadP3 


LoadP1 


Fig. 5. Diagram of simulation grid 


Known parameters to simulate include: source and transformer in Table 1, transmission lines in Table 2, electric load at bused in 


Table 3 and PVs in Table 4. 


Table 1. Parameters of source and transformer 


Type Parameters 

Grid Rated voltage: 110 kV; Short-circuit power: 5000 MVA; Reactance/Resistance: 
Power Voltage ratio: 110/35 kV; Rated power: 25 MVA; Impedance: Z=10%; 
transformer Reactance/Resistance: 20 


Table 2. Parameters of transmission lines 


Name _ Sectional Type Length | Name Sectional Type Length (km) 
area (km) area (mm7?) 
(mm_”) 
P1-P2 183 Pirelli-twisted 19 12 P3-P4 77.3 Pirelli- 8 
strands twisted 7 
P2-P3 111 Pirelli-twisted 7 20 P2-H1 111 strands 30 
P4-H1 49.5 strands 11 H1-H2 49.5 6 
P4-P5 34.4 10 
Table 3. Parameters of electric load at buses 
Name Apparent cos@ Type Name Apparent cos@ 
power power 
(MVA) (MVA) 
LoadP1 2 0.85 80% LoadP5 2 0.85 80% constant 
LoadP2 1.5 0.8 constant LoadH1 3.6 0.8 kVA, 20% 
LoadP3 1.2 0.8 kVA, 20% LoadH1 15 0.8 constant Z 
LoadP4 2.5 0.8 constant Z 
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Table 4. Parameters of PVs 


Type Parameters 

A Sharp NE-Q5E2U Pmppst=165 W; Vmppste=34,6 V; Imppste=4,77 A; Vocst=43,1 V; Iscste=5,46 A 

panel 

PVs Power at standard test condition: 3960 kW; An array: 80 panels in series; 300 
arrays in parallel 


Simulation results representing some operating states of PVs to evaluate bus voltage and power flows in whole grid are shown 
in Fig. 6, Fig. 7 and Fig. 8. 
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Fig. 6. Without solar irradiance (G=0) 
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Fig. 7. Power of solar irradiance is 500 W/m? 
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Fig. 8. Power of solar irradiance is 1000 W/m? 


From simulation results in Fig. 6, Fig. 7 and Fig. 8, the change of power flows in transmission lines is represented in Table. 5. 


Table 5. The change of power flows in transmission lines 


Power flows G=0 W/m? G=500 W/m? G=1000 W/m? 
(kW+jkV Ar) 
Transmission lines 
P1-P2 10520+7906j 8575+7644j 6688+7461) 
P2-P3 4772+3294j 3616+3191; 2484+3116j 
P4-H1 -39-29} 653-89} 1344-117; 
P4-P5 1679+1005; 1689+1010j 1698+1015; 
P3-P4 3650+2435j 2546+2399j 1450+2371; 
P2-H1 424043169} 3523+3151j 2822+3150j 
H1-H2 1173+861j 1178+864; 1183+868j 
PO-P1 12295+9914j 10341+9454j 8448+9115; 
Total power loss 881,8 +1973,3; 668,5+1295,2} 519,7+920,6j 


Results in Table 5 showed that power flows in whole grid change very much although designed power of PVs at standard test 
condition only has 27.7% in total load’s possession with the participation of PVs only at bus P4. With the integration of PVs at the 
end the transmission line, power mobilized from source and power transmitted through the transmission line near source 
decreased very much. When G increased from 0 to 1000 W/m2, the integration of PVs helped to decrease total power loss in whole 
grid upto 51.1% (from 2161.4 kVA to 1057.2 kVA). Moreover, power flows going through P4-H1 transmission line changed: 
direction of active power reversed (increased 3446%) and reactive power increased 403%. From above comment, it must be noted 
when choosing sectional area for P4-H1 line to adapt requirements of transmitted power with the participation of PVs at bus P4. 

Working the same as above operating states in Fig. 6, Fig. 7 and Fig. 8, simulation process will continue with other states 
corresponding to T=30°C and 200 W/m? for each step of G. The change of voltage at bus H1, H2, P2, P3, P4 and P5 is represented in 
Fig. 9. 
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Fig. 9. Change of bus voltage (T=30°C and G varies) 


Simulation results in Fig. 9 showed that the integration of PVs helped to improve voltage quality of all load buses whenever 
having solar irradiance. When G=0, values of voltage at all load buses were below 95%Urated. When G increased and reached to 1000 
W/m?, values of voltage at almost load buses noticeably increased higher than 95%Uratea (except P5 and H2 buses are near 
95% Unrated). It means that PVs can help to improve voltage quality for whole grid and it needs to have participation of PVs at more 


buses to continue higher improvement. 


CONCLUSION 

The important contribution of this paper is to show the effect of PVs to operating parameters in distributed grids, including power 
flows in transmission lines and bus voltage. This paper showed the power structure and detailed control functions of a PVs. 
Moreover, characteristic of PVs was represented to provide some evaluations to dispatchers about the integration of this generator 
to operating parameters. These contributions can help to give out solutions to hold power balance and improve voltage quality for 
whole grid. 

Carried out by ETAP software, power flows in transmission lines and voltage at all buses at any operating state can be 
determined exactly. By using simulation tool, all working states of loads and variations of G can be evaluated in detail to make 
operating plans or forecast the development of the grid clearly. 

The integration of PVs can be continued to evaluated when having participation of other generations such as small hydro plant, 
wind generator, capacitor. Operating parameters of electrical grids can be improved if these distributed generations are 
implemented. It means that optimal locations and optimal capacity for each generation are next research directions to have the best 


effective operation. 
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